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Abstract 

Random Ising systems on a general hierarchical lattice with both, random 
fields and random bonds, are considered. Rigorous inequalities between eigen- 
values of the Jacobian renormalization matrix at the pure fixed point are ob- 
tained. These inequalities lead to upper bounds on the crossover exponents 
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Despite the many years of research and the large number of researchers working on 
the subject, the study of critical behavior of random systems has led to only few exact 
results. On the other hand, some of these results |1],|2[ played an important role, especially 
in the context of the random field problem. In a recent study we considered a random 
bond Ising system on a general hierarchical lattice, where the renormalization group (RG) 
transformation is exact and obtained inequalities concerning the eigenvalues {A;} of 
the Jacobian renormalization matrix, at the pure fixed point. The purpose of the present 
study is to show that similar inequalities can be obtained if random fields are included. In 
contrast to the case of random bonds and zero fields, correlations are now generated by 
the renormalization flow. Nevertheless, it appears that these correlations are, first, confined 
to the fields so that the distribution of bonds is left uncorrelated, and second, restricted 
to nearest-neighbor (nn) correlations. It is important to emphasize that these short-ranged 
field correlations are generated by the RG transformation even if one assumes no correlations 
to begin with, and that the range of correlations does not increase under the transformation. 
Our results are relevant to real lattices, since some approximate RG schemes on real lattices 
are in fact exact RG schemes on hierarchical lattices (Migdal-Kadanoff 0,^] and others 
and since it is believed that the critical behavior of an Ising system on a real lattice can be 
mimicked by that behavior on a properly chosen hierarchical lattice [p|fl |lO|,|lT| . 

We consider a general hierarchical lattice described schematically in Fig. [I]. The shaded 
area shown in (a) consists of a set of lattice points where some of the pairs are joined. In (b), 
a typical shaded area is represented. The solid lines are bonds to be iterated in constructing 
the lattice while the dashed ones are not to be iterated. The bold lines represent the 
possibility for some bonds to be strengthened, namely, multiplied by some constant. The 
random Ising system is represented by the dimensionless Hamiltonian 

—H = JijViVj +^Zh i a i) (1) 

where refers to connected pairs only. All three types of bonds of Fig. |](a) then carry a 
coupling Ji| (for the bonds joining sites a and 0), while each site carries a field h 1 ^ . (Note 
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that one of the members of the pair af3 may be either 1 or 2.) 

The renormalized couplings and renormalized fields are given by 

Jv> = fj{J3,tg] (2a) 

and 

hi = h i + '£f H {J3 h h£} 2 = 1,2 (2b) 

respectively, while Zi is the coordination number of the site % on the renormalized lattice. 
Both, fj and fh, depend only on couplings and fields within the rescaling volume associated 
with the pair of sites (1,2) (the shaded area in Fig. 0(a)). Eq. (|2a| ) implies that Jy and J\ m 
are not correlated if the pairs (i,j) and (l,m) are not identical. This does not hold for the 
renormalized fields. Due to the sum in Eq. (^FJ) over nn sites on the renormalized lattice, 
it is clear that even if there are no correlations to begin with, correlations are generated by 
the RG transformation, between fields on nn sites and fields and couplings on a site and a 
bond attached to it. (For example, in Fig. 0(b), the following pairs are correlated: (hi, /t 2 ), 
( J12, hi) and (J12, h 2 ).) It is easier to deal with such correlations by considering a bond-field 
Hamiltonian of the form 

-H = Yj [ J ij a i a j + hijfa + (T j)] ) (3) 
(id) 

in which the random variables are the couplings Jy and the bond-fields hy. The set of RG 
transformation equations is now given by 

J12 = fj{Jl% h^p} 

and 

Eqs. (^), therefore imply that non of the two couplings jy or hij, is correlated with any of 
the two couplings J\ m or hi m , if the pairs (i,j) and (l,m) are not identical. 
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(4a) 



In terms of the bond-fields, the site-fields are given by 

hi=l^2hij. (5) 

z m 

where indicates that the sum is over all sites j connected to i. A similar bond-field 
Hamiltonian (§) was already used in the past |]l2, r3|,14j| , only with the additional term 



h\j((Ji — Note that it is necessary to include the dagger fields only if one assumes 
that the site-fields are initially uncorrelated. Since, however, the initial state of non corre- 
lated site-fields is not preserved by RG transformation and nn correlations between site-fields 
are developed, there is no reason to start with uncorrelated fields on the sites. 

We assume the existence of a ferromagnetic fixed point at {J a p} = J* and {h a p} = 0. 
We denote the departure of J a p from J* by 5J a p and define the moments as 

r« = ((5J aP Y(h aP y). (6) 

Clearly at the fixed point r*- = 0. We will be interested in the recursion relations of the 
moments near the pure fixed point, 

Tij = Gij{Ti m }. (7) 
The recursion relations above are obtained from the recursion relations for the local couplings 



= E fir) * &Jop + \y, (jw) (^) + \ E (fr ) * ^ + • • • (8a) 

and 

(9fh\* ( d 2 fh \* 

^ V + E ox L $JapKp + . . . , (8b) 

where (...)* denotes evaluation at the pure fixed point. Note that although we are interested 
only in the expansion of to first order in Ti m , we still need, in principle, orders higher 
than linear in Eqs. (|^) above. Note also the terms missing in the expansions due to the 
different parities of J and h. The renormalized moments are given by 
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ij 

(«,/3) 



\fdfjY' 




\( dfh VI 


\dJaf3j _ 







(q,/3) Z,m 



(9) 



where clearly, in the last sum on the right end side of the above, I + m > % + j. Also the 
parity of m in the sum must equal the parity of j. The Ay^s with / + m > i + j always 
involve derivatives higher than the first of at least one of the /'s. We arrange next the T^'s 
using a single index 



(10) 



with 



k = ,y/v o +3 + 1. 



This brings Eq. (^) into the standard matrix notation 



(11) 



(12) 



It is not difficult to show that if k\ corresponds to (i, j) and k 2 to (/, m) then / + m > i + j 
implies hi > k\. This means that the matrix A is block-triangular (Fig. |2|). Consider next 
one of the blocks along the diagonal of A. From the expansions (|8|) and @ it follows directly 
that the only contribution to (5jy(h)i in (5J) l (h) m such that / + m = i + j, is the one with 
/ = i and m = j. The final conclusion thus is that the matrix A is triangular, so that 
its eigenvalues are just its diagonal terms. The eigenvalues of the Jacobian transformation 
matrix are thus 



A 



\fdfjY' 




\( dfh Y] 


\dJap) 







(13) 



This leads now to a number of interesting inequalities: 
(a) All eigenvalues are positive, 



Xij > 0. 



(14a) 
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(b) All eigenvalues are ordered, 



(14b) 



(c) All eigenvalues obey a convexity condition, 

Kj^kj > Aj+fcj and XijXik > \.,j+k- (14c) 

(d) All eigenvalues obey 

(Aij) 2 < \i + kj+i\i-k,j-i, (14d) 

where = — i, . . . , i and I = —j, . . . , j. 

Proof: In a recent paper |J we have considered random bond Ising systems for which 
the subset {Ajo} is considered. There, we have already proven properties (a)-(c) and our 
proof here will follow the same line. 

Properties (a) and (c) are proven by showing that 

dfj 
dJ a p 

and 

dfh 

dh a/3 

We have to consider then the specific transformations generated by 



(J*,0)>0 (15a) 



J*,0)>0. (15b) 



-7Y = lntrV w , (16) 

where tr' represents trace only over the subset of spins {<J a } internal to the rescaling volume, 
not including the external spins o\ and o-i- The renormalized couplings and fields given by 
Eqs. (PS), can now be written in the forms 
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J\2 — — ^tri2 



cr^H , (17a) 



hu — — ^tr 12 



(d 1 + a 2 )H , (17b) 



where tr i2 indicates trace over the two external spins o\ and <r 2 . The derivatives of Jyi with 
respect to J a/ 3 and hu with respect to h a p are thus given by 

dJ 1 

— — (J*,0) = 7tr 12 ((Ji(72)((J a (7^)i2 (18a) 

CJ a/3 4 

and 

<9/i 1 

— -(J*, 0) = -trxa^! + <7 2 )<<7 Q + ^> 12 , (18b) 

where (. . .) 12 is the average with respect to Ti with a x and cr 2 held fixed. In calculating the 
above derivatives at the pure fixed point, we use the following symmetry properties of the 
system, 

(o"aO"/3)+ + = {a a a p )*__, (19a) 
(oaO"/3)+_ = (cr a vp)*_ + , (19b) 

= (i9c) 
M* = -M* (i9d) 



to obtain 



d J 12 1 , 



-(J*,0) = - (<7 Q ^);_] (20a) 



and 



|ii (J -,o) = i[^ + ^; + ]. (m) 

Here the sign indices specifically indicate the state of the spins <j\ and <r 2 and the * indicates 
that the average is with respect to the pure fixed point Hamiltonian 
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-w = rY t wy (2i) 

Now, according to the GKS inequalities [15,16], if all the many-spin couplings J a = 
h a , J a /3, ... in a general Ising system are positive, all the many-spin correlations (a a) = 



(a a ), (cr a (j/3), . . . must obey (a a) > 0. Using Eqs. (p20|) , the averages are taken with respect 



to the pure ferromagnetic Hamiltonian (21), where the two external spins of each of the 
rescaling volumes, which are held fixed, serve effectively as local fields. When these effective 
fields are held both positive, the GKS inequalities hold, so that 

M* ++ >0 and (cr a (Tp)l + > 0, (22) 



which is enough already to prove Ineq. ( 15b ). Ineq. ( |15a| ) can be easily shown to hold 



using, in addition to the GKS inequalities, other rigorous inequalities, just recently proven 
|T7| , also concerning the many-spin correlations in general Ising systems. It states that if 
all the many-spin couplings J a are positive again, the absolute value of all the many-spin 
correlations {a a) does not increase when the value of any of the couplings is reduced, taking 
any value in the interval [—J a, J a}- According to this, it is clear that under reversal of the 
+1 state of any of the two external spins, the many-spin correlations can not increase. So, 
we arrive at the conclusion that 

(cr a crp)* ++ > (a a ap)\_ (23) 

and Ineq. ( |15a| ) is also proven. This completes the proof of properties (a) and (c). 
We turn now to property (b). here we need to show that 



dfj 
dJ a p 

and 

df h 



(J*,0)<1 (24a) 



(J*,0)<1. (24b) 



dh aj3 

at any finite temperature. But, referring to Eqs. (|20|) again, it is clear that 
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^[((Ta^p)*++ - {<Ja<?p)*+-\ < ^[\((? a epY ++ \ 



ik^>;_|]<i 



(25a) 



and that 



+ <^>; + ] < \[\KY ++ \ + < 1, 



(25b) 



while the equality sign can hold only at zero temperature. This proves property (b). 

We are left now with property (d). Here we use the more general definition of a scalar 
product, (u, v) = J2i w i u i v i where Wi > and the corresponding Schwartz inequality, 
which reads (J2i WiU*v «) 2 < \ui\ 2 J2j w j \ v j\ 2 (here is the only place where the * repre- 

sents complex conjugate). We replace, next, the sum over the single index i with the double 
index (af3) and identify 

"/ox \*~l r r/a-f \ *~1 s 

(26a) 



\fdfjY' 


r 


\( dfh \] 


\dJap) _ 




\dKf}) _ 



\fdfjY- 


p 


\( dfh Y\ 


\dJ a p) _ 




\dh afS ) _ 



U af 3 = 

with r = 0, . . . , i and s — 0, . . . ,j, 

V af 3 = 

with p — 0, . . . , i — r and q — 0, . . . , j — I, and 
\<> obtain 



(26b) 



\fdfjY~ 


i—r—p 


\( dfh X] 


\dJ al3 ) _ 




\dh a/3 ) _ 



3-s-q 



(26c) 



i-r+p,j~s+qi 



(27) 



where we have used the fact that all partial derivatives are real and positive. All is left now 
is to denote r — p = k with k = —i, . . . , % and, similarly, s — q = I with / = — j, . . . , j, which 
completes the proof of property (d). 

In addition to that, denoting by mj < 1 and m h < 1, the maximal values of dJij/d.J a p 
and dhijjdha.fi respectively, we obtain 



\j < himj k m> h \ with 



k — 0, . . . , % 
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and I — 0, ... ,j, 



(28) 



so that we have also proven that the number of relevant interactions at the pure fixed point 
is finite. The only case of which the equality sign hold is the diamond hierarchical lattice 
(DHL) P|JIB|], where all bonds are equivalent. From Ineq. (pBD follows an inequality for the 
crossover exponents: 

9ij < H ;— : <1 i+j = 2,3, (29) 

m Aio 

where 0^ = yij/yio, yij = In Ajj/ln5 and b is the rescaling factor. The condition for criticality 
of the pure fixed point is max(A2o, An, A02) < 1, while else, we expect a random critical 
point with a different set of critical exponents. It is interesting to note that it was just 
recently shown |19[ that even for the random bond Ising system, the Harris criterion for 
pure criticality p0[] , a p < (a p being the specific heat exponent) is equivalent to the 



obvious requirement, A20 < 1, only in the special case of the DHL. In the more general case, 
it was shown that a p < 4>20 so that the Harris criterion is only a necessary condition for pure 
criticality to hold and counter examples where a p < and 020 > have been presented. The 
analogous result for the random system is ■y p < O 2 (jp being the susceptibility exponent of 
the pure system) but since 7 P turns out always to be positive, the random field is always 
relevant at the pure critical point. 

We wish to conclude by emphasizing that the inequalities proven here hold not only for 
exact RG transformations on HL's but also for all other renormalization schemes (such as 
MK scheme |5]||) in which the renormalized couplings or fields are not correlated or even 



in cases where it is clear that the correlations are not important [21 
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FIGURES 




(a) (b) 

FIG. 1. A general hierarchical lattice is described schematically. In (a), the shaded area consists 
of a set of lattice points, a, (3, . . ., where some of the pairs are joined. In (b), a typical shaded area 
is represented. The solid lines are bonds to be iterated in constructing the lattice, while the 
dashed ones are not to be iterated. The bold lines represent the possibility for some bonds to be 
strengthened, multiplied by some constant. 
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FIG. 2. A general schematic description of a "block-triangular" matrix is shown. Bi, . . . , B4 
represent the blocks along the main diagonal. The gray area marked with X indicates the presence 
of non zero matrix elements while in the area marked with 0, all elements are zero. 
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